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Abstract 



We build a simple cosmological model by means of a parabolic 
Lemaitre-Tolman-Bondi (LTB) metric up to a redshift z ~ 0.4, an hy- 
perbolic Friedmann metric for z ~ 0.4 up to the scale where dimming 
C^) \ galaxies are observed (z ~ 1.4) and a bulk spatially flat metric up to 

• the last scattering surface. Following Wiltshire, by taking into account 

the different rate of clocks for an observer at the centre of a parabolic 
■ LTB spacetime with respect to a one in the hyperbolic Friedmann met- 

Q> | ric, an "apparent" negative deceleration parameter is perceived by the 

observer at the centre of LTB, provided that all the regularity con- 
ditions are imposed and the past null sections of the LTB and the 
hyperbolic Friedmann metrics are identified. As a result, a first order 
Hubble law emerges at low redshifts. A parameter K arises driving 
the deceleration parameter perceived by the central LTB observer. Fi- 
nally, we obtain that a negative value for the deceleration parameter 
is compatible with the observed energy-density at our present epoch. 

Keywords: Dark energy, deceleration parameter, Lemaitre-Tolman-Bondi 
metric, clock effects 

PACS Numbers: 98.80.-k,98.80.Jk,95.36.+x,04.20.-q 

1 Introduction 

Supernovae type la (SNIa) observations of the past decade seem to indi- 
cate an accelerating universe ([HE]). In the standard approach with the 
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Friedmann-Lemaitre models (FRLW), an accelerating universe invokes the 
presence of a large amount of the so called dark energy. In the FRLW 
picture, this dark energy is given by the cosmological constant. The dark 
energy represents a puzzle and perhaps the biggest problem in the mod- 
ern cosmology. In fact, a direct detection of a cosmological constant is still 
lacking. In the last decade, many attempts have been made (see [31 [19] 
and references therein) to obtain physically sensible models predicting a 
negative value for the deceleration parameter. Some authors (see for exam- 
ple P3J HU HS1 Q21 HH1 [H]) showed that inhomogeneities can generate an 
accelerating universe by using LTB metrics (see [201 ED 122]), but several 
conditions must be imposed (see O [17] ) in order to build regular physically 
viable models. In particular, in [16] it is shown that LTB metrics can mimic 
the distance-redshift relation of the FRLW models at least at the third order 
in a series expansion with respect to the redshift near the centre where the 
observer is located. More generally, in the LTB solutions, apparent accelera- 
tion in the redshift-distance relation seen by a central observer can be shown 
to coexist with a volume average deceleration on a spacelike hypersurface 
(see [23] ) An accelerating universe can also be builded by averaging inhomo- 
geneities (see [3 SI El [71 O HOl HH H2l HU ) by means of the techniques 
depicted in El [7] . For a review on inhomogeneous cosmological models see 
[24\ 125]. Particularly interesting is the idea developed in Wiltshire's papers 
9. 12J. In these papers the dimming of the distant galaxies is interpreted 
as a "mirage" effect. This effect is due to the different rate of cloks located 
in averaged not expanding galaxies, where the metric is spatially flat, with 
respect to clocks in voids where the spatial curvature is negative. To a neg- 
ative spatial curvature can be associated a positive quasilocal energy. This 
gravitational energy is not local, according to the strong equivalence princi- 
ple. This approach seems to be very promising. 

In this paper we adopt the point of view present in [9l Q2]. In particular, 
we are interested in models that take into account the observed inhomo- 
geneity of the universe at least up to a redshift z ~ 0.1 by means of LTB 
metrics. We apply the reasonings of [11], that, as claimed by the author, 
represent the first level step in solving the fitting problem (see |26j). We use 
for our local near universe a parabolic LTB metric, instead of a "dressed" 
(see [7]) hyperbolic Friedmann metric. In particular, in the spirit of [11] . 
we obtain a simple (crude!) model. Our starting point is the consideration 
that dimming galaxies are located at hight redshift (z ~ 0.5 — 1.4) [27] . 
Therefore, we build a model where, up to z ~ 0.4, the universe is filled with 
a parabolic LTB metric, for z > 0.4 with an hyperbolic Friedmann met- 
ric and for z > 1.4 with a bulk spatially flat Friedmann metric, according 
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to WMAP satellite data [28]. Similarly to [TT], we neglect the coupling of 
the dynamics of the different scales depicted above. In this way, the SNIa 
dimming galaxies are in a region where the clocks are ticking slowly with 
respect to the ones of an observer located at the centre of a parabolic LTB 
metric. Within this crude model, firstly we obtain (see a linear Hubble 
law for small redshifts. Further, we show that a negative central decelera- 
tion parameter arises, provided that the necessary boundary and regularity 
conditions are imposed. Furthermore, we show that a negative central de- 
celeration parameter emerges only for values of the density compatible with 
the ones actually observed (underdensity) . Following Wiltshire, since the 
SNIa observations are made along the past light cone, we identify the null 
sections (9, (f) = const.) of the parabolic and the hyperbolic "dimming" zone. 
It should be noticed that cosmological models with only two metrics can be 
exhaustively found in |30j . 

In section 2 we present the metrics composing the model together with the 
initial an regularity conditions. In section 3 we present our model. In sec- 
tion 4 we obtain the linear Hubble law for low redshifts. In section 5 regular 
solutions are discussed. In section 6 we obtain a distance-redshift relation 
together with the expression for the central deceleration parameter. Sec- 
tion 7 collects some final remarks and conclusions. Finally, the appendix is 
devoted to the study of the matching conditions. 



2 Initial and regularity conditions 

The starting point of our simple model is the consideration that the dimming 
galaxies have been found at hight redshifts {z ~ 0.5 — 1.4). Consequently, we 
assume that the universe with z ~ 0.5 up to some units, is represented by an 
hyperbolic Friedmann metric with negative spatial curvature. In appropriate 
coordinates, the metric can be put in the form 

ds 2 F = -dt 2 + a(t) 2 [drj 1 F + sinh 2 ^ d£l 2 ], (1) 

Hit = ^_ 3 (sinh£ - £) , 

2(1 - 

a(t) = a% — - (coshf — 1) , 

a(t) 2 H 2 (l -U) = 1, 

where is an initial density parameter, Hi an initial Hubble constant and 
Hi an initial expansion factor to be specified. An observer in the portion of 
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universe given by dU) measures the observables by means of the comoving 
time t. With respect to this time, an observer in ([1]) measures an Hubble 
flow with a time dependent Hubble constant H given by 

adt ^(cosh£- iy ' y ' 

In what follows we adopt the simplifying assumption that the dynamics (see 
of the pieces composing our model are independent. This is a "crude" 
assumption, but nevertheless it gives the possibility to study the role of 
inhomogeneities and structures in the distance-redshift relation. According 
to WMAP data [28] . over the "dimming" zone, the universe can be modelled 
with a bulk Friedmann metric with zero spatial curvature, i.e. 

ds B = -dt\ + a 2 B (t B )(dr) 2 B + rf B dU 2 ), 

2 

ob(*b) = a,Bi f^rj ' ( 3 ) 

where in $5$) we have assumed a dust model. With ([3]), the Hubble flow H B 
is given by 

H B = = — . (4) 

as dts 3ts 

Finally, we assume to live at the centre of a parabolic (vanishing spatial 
curvature) LTB spacetime up to a maximum of z ~ 0.4. The metric is 

ds 2 TB = -dT 2 + R 2 drj 2 + R 2 dQ 2 , (5) 



where, as usual: 



47TP= R^> (6) 



mv) = (^y [T - Yom i, (7) 

and subindices with comma denote partial derivative. The arbitrary function 
Yq(t]) is often called "bang function" and is usually interpreted as a big-bang 
singularity surface. The arbitrary function M{rj) represents the gravitational 
mass inside a volume of radius rj and p the local density. 
Regularity conditions must be imposed to ©-([7]). First of all, the density 
([!]) must be positive everywhere and finite at the centre 77 = where the 
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observer is placed and R(T, n) must be vanishing at r] = 0, V T. 
Mathematically (see [291 13"T] : 

p(T,rj^0) = finite, (8) 
R(T,n^0)~nf(t) , R v >0 Vr?,i, (9) 
R(r], T) > V T and n > , M(r? -> 0) ~ n 3 . (10) 

Further, no trapped shell singularities must arise in the TB zone, i.e. 

n > 2GM(n). (11) 

The absence of trapped shells is automatically satisfied near the centre, 
provided that conditions (f9j)- ([T0j) are imposed. We also impose (see [17] ) 
the condition 

Yo, v (v = 0) = 0, (12) 

to avoid potential problems at the centre, such as the "weak singularity" 
discussed in [TTJ, although the severity of this problem is debated (see [32], 
[33]). In section 6 we show that the presence of this "weak singularity" does 
not affect the central deceleration parameter. 

We consider now the initial conditions that we must impose to the metrics 
(TjQ) , ([3]) and ([5]) . If we consider (see [H] ) an early time such that Qi is close 
to unity, all the three scales depicted above must be matched at that early 
time. As a result Hi ~ a^i- 

The same condition for ([7]) it gives R(T{, n) = nHi. This condition, in terms 
of the bang function Yq reads: 

To study the condition ([T2I) is more useful to write Yo(n) as follows: 

«--tA' (14) 

In this section, to represent the dimming "zone", we used the hyperbolic 
Friedmann solution. This choice leads to simple computations, while an 
hyperpolic LTB spacetime generally does not allow to simple analytic ex- 
pressions. For example, the volume average (|15p for unbound LTB metrics 
is generally not at our disposal in an explicit workable form (see [3]). Fur- 
ther, note that the approximation used in this paper can be justified both 
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physically and mathematically. In the appendix we show that the continu- 
ity of the first fundamental form can be achivied on a non-comoving thin 
shell, ensuring that the spacetime is connected. Further, the Stephani met- 
ric (see [31]) could be used as a comoving thick shell located between the 
two transition zones composing our model (see the appendix). It is worth 
to be noticed that to calculate the distance-redshift relation only the con- 
formal null sections of the metric come in action, since the astrophysical 
observations are performed along the past null cone. 



3 The model 



The further step in our study is to write the relations between the differ- 
ent scales composing the model. To this purpose, it is observed a broadly 
uniform Hubble law (see [12J and references therein). Therefore, following 
Wiltshire, we impose the equality of the Hubble flow of the three scales. 
Concerning the metric ([5]), we can calculate, in the spirit of the Buchert 
scheme [5], a volume average expansion up to some scale r/D, where dim- 
ming galaxies come in action. For the proper volume Vd and the expansion 
9, we read: 



'ID 



V d = Att R r ,R z dr], 



<V>D 



a D (T) 



d 



(lnVb), M D = M( VD ) 



dT 

(T-Y (r, D )) 



(15) 
(16) 
(17) 



XTi-Yo(r, D )) 

As a result, for the averaged Hubble flow Ttb i n the LTB sector, we have 

< 6 > D 2 



Htb 3 3[T-Y (r, D )Y 
After equating ©, @ and ([15]) . we get: 
t B = T 



T(0 



Y (vd) 



(cosh£ — 1)" 



3ffi(l 



sinh^ 



+ Y (rj D ) 



(18) 

(19) 
(20) 



Apart from the constant in the right hand side, the expression (|20p for the 
time delay is the one found in [TT]. This is not a surprise because, after 
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averaging a parabolic LTB metric, we obtain an Hubble flow, equation (|18p . 
that, apart from the constant Yq(i]d), is the same of the bulk spatially flat 
metric. The constant Yq{t}d) takes into account the cut-off made in our 
model. Further, the term Yo(wd) appears as a translational factor and thus 
does not enter in our analysis. We only mention the fact that such a constant 
can give a correction (positive if Yq(j]d) > 0) to the age of the universe. For 
the lapse function J(£) = we get: 

i(t\ - 3 (l + coshg) 

J( ^- 2 (2+ cosier (21) 

Formulas f)20[> - (|2 1 1) describe the different rate of clocks between the parabolic 
TB observer and an hypothetical observer placed where the galaxies are 
dimming. Obviously, for the reasonings above, formula (|2ip is exactly the 
one found in [Tl], but is expressed in a different background. 
If we want to describe the dimming of distant galaxies, we must to relate 
the metrics (pQ) and ([5]) on the past null cone, where the SNIa observations 
are performed. To this purpose, the radial null sections (6, (ft = const.) of 
(PQ) and ([5]) must be the same, i.e. ds 2 F = J 2 ds TB . As a result, the following 
equations hold on inward radial null geodesies: 

- dT = -dm, (22) 
a 

R v dn = -dT. (23) 
Therefore, our related metric in the LTB inhomogeneous spacetime is 

a 2 

dsl B = -dT 2 + -^dnl + R 2 dQ 2 . (24) 
J z 

For a more complete discussion regarding the matching conditions see the 
appendix. 

For the metric (|24|) , we can define a radial observed (by the central observer) 
Hubble flow with 

H ob = 1 -^ = 4,(^)- (25) 



a dT dT \ J 

It is in terms of (|25p that we measure the Hubble flow. Concerning the 
luminosity-distance d^z), for the metric (f24l) we get (see [T5 | l35 | l36 l 137 } 138 } 
[39]): 

d L = B (l + z) 2 ,B 2 = ^ = R 2 , (26) 

as io 

where Slo is the solid angle subtended by a bundle of null geodesies diverging 
from the observer and Sq is the cross-sectional area of the bundle. We must 
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to integrate the relevant field equations for our purpose. First of all, by 
means of ([I]), ([2T|) we can integrate equation ([22]) along the past null cone. 
We get: 

VF = £o - £ , £ < £o, (27) 

where the subscript "0" denotes the actual time related to the central TB 
observer. Following Celerier [15], for the metric (|24p we can express the 
observed redshift z in terms of the time parameter £, i.e.: 



1 



1 



(1 + z 



_d_ 

dT 



!)' 



(28) 



l + z 



(29) 



and therefore, by means of (|22h and integrating backward starting from 
z = 0, we read: 

J a 

Since the expression (|2"9~|) has been evaluated along the past null cone, it 
is the same found in [11]. Nevertheless, since inhomogeneities have been 
taken into account in our model, we expect corrections with respect to the 
picture of the paper [11], in particular in the relation distance-redshift di,(z). 
In evaluating the functions entering in the relation di{z), we need of the 
difference T — T, with T the actual time. The equation ([29]) permit us 
to express cosh£ in terms of the redshift z. As a result, after noticing that 
J7o = i +C Q S h g an d with the help of ([1]) and (|20p . we obtain: 



(cosh£ 



1)1 



A 



T = A 

(1 + cosh£o) t 

n (2 + n 2 ) i 
H {i-n f 2 (2 + ^ ) 2 



(cosh £ - 1)2 
(1 + cosh£)2" 



(30) 



where Hq is the measured Hubble constant given by (|25[) and calculated at 
the present time £o (or To), i.e. 



/To — 3T/o 



2 + ^q 

(2+n ) 2 . 



(31) 



Finally, along the past null geodesies we have = My^, and thus, 
thanks to (f23|) . expression ([6|) becomes: 



P(T) 



' 4vri?(T) 2 ' 



(32) 



S 



Obviously, regularity of ([32]) for T — ► To requires that: 

M(T^T ) ~ (T -T) 3 + O (l), (33) 
R(T^T ) ~ (T -T) + O (l). 

4 Zeroth order solution: the linear Hubble law 



To complete our model, we must integrate the equation (|23|) . Obviously, 
because of the partial derivative of R in the left hand side, this equation 
cannot be integrated in this form. Nevertheless, the equation f|23|) can be 
easily integrated as follows. Firstly, we write: 

f dR(T, V (T)Y 



drj J \ dT 

and the equation f)23[) becomes 



dT = dR = R >v dq + R, T dT, (34) 



dT = W^Ty (35) 



Further, from ([7|) we obtain R % t = \ Jr , and as a result 



dT = _ g _ , (36) 

2GM _ i 
R 1 

We will explain the general strategy to integrate the equation (|36p in the next 
section. First of all we are interested in the first order calculation of di{z). 
The first condition to impose is R > 2GM, that is equivalent, near the 
centre, to the first of conditions (jlip and ruls out trapped shell singularities. 
It is worth to noticing that the conditions (|33|) near the centre are sufficient 
to satisfy the condition above mentioned. The zeroth order approximation 
for (|36h emerges when R >> 2GM. This extreme approximation it gives 
the correct first order of d^(^). Consequently, after integrating with the 
appropriate boundary condition (R(Tq) = 0), we have: 

R{T,T](T))=R{T) = T -T + o(l). (37) 

With the help of (j29]l . (|30]h we read: 

d L {z) = ^- + o{z), (38) 

that is the well known linear Hubble law for low redshifts. Therefore, the 
zeroth order of our model it gives the observed distance-redshift relation for 
z « 1. 
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5 Exact regular solutions 

We can rewrite equation ([36]) as: 




mr) = -l ( i a/^ )dT. (39) 



To integrate ([39]) . we write J 2< ^ M = F(T), being -F(T) a regular differen- 
tiable function. With the conditions (|33p we must impose: 

F(T) G (0, 1) , fl(T) > 0. (40) 

In this way, after fixing an ansatz for F(T), we can integrate the equation 



39]) . obtaining: 

T 



R(T) = T - T + f F(T)dT. (41) 
After solving the equation (|4ip for R(T), M(T) is given by 

M(T) = ^F{Tf. (42) 



Further, we impose the condition (|12f) . that in terms of equation (|14h seen 
as a function of T becomes: 

io,r(To) = 0. (43) 

Finally, the relation between rj and T along the past null cone is obtained 
by inverting the equation ([7|) with the help of (fl~3|) i.e. 

2 



/9GM\3 



Rhl— T + Tj 

w 9GM 



(44) 



In the following (in particular for the distance-redshift relation (|48|) ) it is 
essential the behaviour of F(T) (and R(T)) near the observer at the centre. 
As a result, any given expression for F(T) must have a taylor expansion 
near the centre fixed by the regularity conditions (|40|) and (|43|) . Hence, we 
can take for F(T) a polynomial expression. Therefore, for our purposes, 
without loss of generality, we can take for T <Tq 

F = ^(T - T) + |f (T - Tf + Q# 3 (T - T) 3 (45) 

« = T o-r-^(T -r) 2 -^(To-r) 3 - 



2K K u y 3i^ 2 
Q Hl(T -T)\ (46) 
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with K and Q adimensional constant. It is worth to be noticed that the 
first two terms in (|45p are fixed by the condition (03]). If we impose that 
the LTB scale is extended up to a maximum of z ~ 0.4, we see that the 
conditions (|40p are satisfied for K > i, with Q of the same order of K. 
Further, note that the limit z ~ 0.4 for the boundary of the TB metric 
only changes the allowed values for K and Q. For example, for z < 0.4 
K > a > i. Furthermore, note that if we impose the condition Yq > for 
T — > To with Yb(To) = 0, we must have: 

K = ~H T , (47) 

while, if we take K < ^HqTq (but positive), then Yq > with Yo(To) 7^ 
0. As we see in the next section, the positivity of Yq, i.e. condition (|4Tj) 
implies a maximum possible value for the central deceleration parameter. It 
should be noticed that if we take a more general expression other than (|45p. 
in order to satisfy all the regularity conditions depicted above, the Taylor 
expansion of F(T) near the centre must be equal to expression (f4"5j) . at least 
for the first two terms. In the next section we show that are exactly these 
terms that enter in the expression for the central deceleration parameter. 
As a final remark for this section, note that in our model, thanks to the 
equation (j42[) . the function F(T) is related to M(T) on the light null cone. 
Further, by means of the equation (|44p we can (at least in principle) find 
M = M(n), i.e. the dependence in terms of n. Consequently, with respect to 
our construction, the FLRW limit can be obtained by setting the particular 
expression for F(T) on the past null cone such that, when expressed in terms 
of n by means of the equation (f4"4"|) . we have M(rf) ~ rf or Yq(j]) = constant. 



6 Central deceleration parameter and observed den- 
sity 

With the help of p9|) . (i30j) we can express M(T),R(T) in terms of the 
measured redshift z. By taking the expressions (|45p . (|46|) and after a Taylor 
expansion near the centre {z = 0), we get: 



+ 



(48) 



4H K(2 + nlf 



8 + 16K 



+ o(z 2 
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The central deceleration parameter go is given by (see [151 El SO] ) 



d 2 

go = ~ H o^2 (M* = °)) + !> ( 49 ) 



that with (1481) becomes 



go 



-3^(5 + 2Qq - 4KSf + 8^o + §KVil + 8-8K 



2K(2 + Q 2 ) 2 



(50) 



If we do not consider the flat Friedmann metric, we could to estimate the 
bulk central deceleration parameter goB at early times by setting £Iq = 1 in 
(f50j) . As a result: 

go* = -\ + ± (5i) 

Note that for K — > 1 it follows that gos — \- 

By taking the asymptotic limit To — > oo (f^o — > 0) we have 

gooo = -1 + ^- (52) 

From (|52p we see that our model is consistent with an accelerating universe 
provided that K > 1. For very large values of K we obtain gooo — — 1- If 
we put in ([50]) an estimate value for Q,q as, for example, < -j^j , we obtain 
that for X — > 1 + again go < 0. More generally, go < at the times for 
which: 

K > _^+f°-f • (53) 

(3O + 40q - 6fi + 8) 

For example, for f^o — ^j, K > 1.0017. Therefore, also for the actual uni- 
verse our model admits a negative central deceleration parameter. By taking 
the limit value for K given by (|47p . (Yo(To) = 0) we obtain a maximum neg- 
ative value for gooo i-e. gooo = — \- As a final step, we consider the local 
density given by (|32p. It is a simple matter to see that, with expressions 
(f45|) . (|46l) . by taking the limit T — > To, we obtain: 

« = P(T„) = ^|_. (54) 

Expression (|54p represents an interesting result: an actual underdensity is in 
agreement with a negative value for go- This is in agreement with the actual 
exstimation predicted by Friedmann models with p c = and y <\. 
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Finally, note that our calculations can be easily changed if we do not impose 
the condition f)12[) . We must only make the following substitutions in the 
equation (f4"5j) : K —> a, K 2 — > b, being (a, b) > (the weak singularity is 
ruled out when a 2 = b). As a result, since the go parameter only involves 
the second order in Tq — T in the expression (|46|) . the central deceleration 
parameter remains unaffected and so also expression f)54|) does not change. 
The so called weak singularity comes in action only at the third order in z 
in the redshift-distance relation. 

7 Conclusions 

Following Wiltshire papers, we build a model for the universe without dark 
energy, by taking into account the observed inhomogeneous universe for low 
redshifts by means of a LTB metric. When the clocks of the observer at 
the centre of a parabolic LTB spacetime are related to the ones placed in 
an hyperbolic Friedmann metric where are located the dimming galaxies, 
an apparent negative value for the central deceleration parameter arises, 
provided that all the regularity and boundary conditions are imposed. We 
attempted to introduce a model in which both a description of the low red- 
shift irregularities and of the large scale homogeneity where present. Since 
a coupling between the dynamics of the different scales composing the uni- 
verse is neglected, the model is a "crude" approximation. Nevertheless, 
our model allows for a negative value of the central deceleration parameter 
that is in agreement with the observed underdensity for the actual universe 
[28] . With respect to the Friedmannian picture of [TT], the use of a LTB 
parabolic metric to describe our nearby universe, permit us the introduction 
of a dimensionless parameter K whose allowed values are in agreement with 
an accelerating universe. In [11] the deceleration parameter runs to zero 
asymptotically from positive values. 

In any case, we have also shown that the so called weak singularity comes 
in action only at the third order in z in the expression ([48ft . and as a result 
the parameter qo does not "feel" such a hypothetical singularity. 
Physically, at the scale where dimming galaxies come in action, we adopt 
an hyperbolic Friedmann metric. In fact, with such a metric is associated a 
positive gravitational energy that is the "source" of the slowing rate of the 
clocks when compareted with the central (parabolic) LTB observer. Accord- 
ing to WMAP satellite data [28] and CMBR [H] over the "dimming" zone, 
we adopt a spatially flat Friedmann metric. As a final remark we cite the 
recent paper [42J where, by means of an investigation of new sperimental 
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data, it seems that we live in a universe with a deceleration parameter near 
to qo ~ 0~, i.e. with a slowing down of the cosmic acceleration. These re- 
sults seem to be in disagreement with the standard LCDM model and could 
encourage the point of view of our work. 

A further development of this paper could be to consider, as in [IQICE]) the 
dynamics by means of the full Buchert [5] formalism by taking into account 
the inhomogeneous structure of the nearby observed universe together with 
the backreaction or to use the full covariant machine given in [6]. Not a 
simple task! 
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APPENDIX 

We now study the matching problem between a parabolic LTB solution and 
an hyperbolic Friedmann solution on a regular comoving surface S. Since 
of the spherical symmetry of both metrics, we can use spherical symmetric 
coordinates £ a on S. Therefore, if we denote with tp^_ the coordinates of the 
Friedmann metric and with tp" the ones of the LTB metric, we have (see 
|35j ) £ Q = ip+ = and, as a result, we can take £° = t = T, £ 2 = 9, £ 3 = <fi. 
In this way the comoving surface S is given by: 

S- = 7] — So = , S + = rjp — So = 0, (55) 

where So denotes the boundary of the LTB metric. The continuity of the first 
fundamental form (see 031111]) on S, i.e. ds 2 s = gapd^d^ (a = (0,2,3)) 
leads to ds 2 _ = ds 2 + , i.e. 

(R)\s = a(sinhr/ F ) |5 . (56) 
For the unit normals n~ and nt we have 

% = R ,T) S l ' K = ^ 5 ]i- ( 5? ) 

For the second fundamental form K a p = (n a .^)| S we have K a p = ^^^fK^ 
The continuity condition on S, i.e. = K^g, becomes 

1 = (coshr/i?)^. (58) 
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From equations (|56p and (|58j) it is evident that is not formally possible to 
match an hyperbolic Friedmann metric with a parabolic LTB one on the 
comoving surface (|55p . 

If we take for the outer metric (+) an unbound comoving LTB line element 

R 2 

ds 2 TB = -df 2 + -g-dn 2 + R 2 dtt 2 , (59) 
f 2 

with f 2 > 1, insted of flSS} we have 

(ij) |s = (.R)| S , T = f+ constant , (/) )s = 1. (60) 

As a result, for the lapse function J we have J = = 1. 
Remember that, since astrophysical observations seem to show a broadly 
uniform "Hubble flow" (see [E]), we must perform the matching according 
to this astrophysical evidence. As a result, the equation (|2ip must be satis- 
fied for the lapse function (J / 1). Hence, as an example, we can think to 
match on a non-comoving thin shell. To do this, we can take for the shell: 

V = Vs(T(t)), VF = n f Mr)), (61) 

where r is the proper time on the shell. Since we have T = T(r), t = t(r), 

the continuity of the pull-back of the metric on the non-comoving shell S it 
gives (on S): 

dT 

Tr = < 62 > 

l-B, (63) 

a sinhr/j s = R, (64) 

£ = J, (65) 



A = Jl + R%f]l, (66) 



B = yjl + a 2 f, 2 s , (67) 

where dot is the derivative of the proper time on the shell (all the expressions 
are calculated on the surface (|6T|) ). The system (|62l) - (i67"l) can be integrated 
in different ways. As an example, from the equation (|62j) we have 




dr= Wl- i? 2 ,,^- (68) 
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Thanks to ([68]) . equation ([63]) becomes: 



The equation ([64]) permit us to express smlir^ in terms of (rj s ,T), and thus 
we can use equation (|69|) to resolve for . Finally, we can integrate equation 
(I68p to obtain the relation r = H(T). 
For the unit normals we have: 

V - = [-rj s R v , R V A, 0, 0] , (70) 
r£ = l-Vfsa, aB, 0, 0] . (71) 

The equations for the continuity of the extrinsic curvature are: 

ij s R tV RR t T + RA = ^f s a 2 a t tsixih 2 r]f s , + 

aB sinhr/j s cosh r]f s , (72) 

rj s R, v f - R tV ri s A + 

r^R^R^T - R, v , v r]sA = 

r)f s ai - aBijfs + rf fs a 2 a^ t - 

Conditions (|72p seem to be incompatible with the equations (|62[) - f)65|) . In 
particular, we have no sufficient number of functions to satisfy the condi- 
tions fl72j) . 

Concerning the matching between the hyperbolic and the parobolic Fried- 
mann solutions, for the first fundamental form we have the equations (162p ~ 
([67|) with R^ — > as in ([66]) , R — > as in ([M]) and J = instead of the equa- 
tion (|65p (with the same J!) and, obviously, with different non-comoving 
surfaces and proper time on the shell. Hence, also for the Friedmann met- 
rics all the reasonings after equation (j67[) are also valid. As a result, at 
least with respect to the intrinsic curvature, the matching conditions can 
be satisfied. It is worth to be noticed that the continuity of the metric on 
the non-comoving thin shell is the minimum requirement ensuring that the 
whole spacetime is connected. 

We briefly discuss another possibility. We can take the spherically simmetric 
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perfect fluid Stephani metric (see PU): 

V 2 

1 

Y{t) d (V 



ds 2 = -D 2 dr 2 + [dr 2 + r 2 dn 2 } 
V = l + ^k(r)r 2 , 



Z?(r,r) = F{t 
fc(r) = Y 



V dr \Y 

2 



C 2 (r) 1 



F 2 {t) 



(73) 



being F, C, Y free functions. The metric (|73[) can be used to model a co- 
moving thick shell (see |45] ) between the LTB and the hyperbolic Friedmann 
zone and so also between the hyperbolic and the parabolic Friedmann met- 
rics. Differently from the thin shell discussed above, the metric (1731) has, at 
least in principle, a sufficient number of arbitrary functions to perform the 
matching by imposing the continuity of the first and the second fundamental 
form, avoiding surface-layer matter, althought in practice this can result not 
easy. However, in this case the calculations of this article could be changed. 
In any case, if this thick shell is "small" with respect to the glued regions, 
one may believe that the deviation from the calculations performed in this 
paper remains "small". In particular, if we denote with 771,772 the boundary 
of the thick shell between the parabolic and the hyperbolic metric, then if 
772 — 771 is "small" , we expect a "small" deviation on the redshift-distance 
relation. Obviously, we could use more general metrics for the thick shell 
than the Stephani one. In a future work we shall consider the model of this 
paper with the introduction of a comoving thick shell. 
As a further remark for this appendix, it should be stressed that, from as- 
trophysical data, emerges the necessity of more scales to describe the whole 
universe we observe. If this is the case, the observed broadly uniform Hub- 
ble flow imposes the equality of the (spatial averaged) Hubble flow of the 
different metrics composing the universe we observe. As shown in this paper 
(see equation (|20j) . the equality of the Hubble flow leads to a clock delay 
effect. If this reasoning is correct, as a consequence, the time flow cannot 
be chosen globally uniform. 

In this paper, the matching is performed by equating the conformally re- 
lated null sections with conformal factor given by J (the lapse function). 
It is worth to be noticed that for the calculation of the distance-redshift 
relation, only the null geodesies come in action. Generally, since only null 
geodesies can probe the cosmological scales under consideration, it seems to 
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be reasonable to impose, at least, the matching between the null sections. 
In this way we have a crude but sufficient approximation to explore the 
role of the observed inhomogeneities on the distance-redshift relation. As a 
final remark, note that, with the introduction of a LTB parabolic metric, 
we have shown that the calculation performed in [llj are compatible with 
an accelerating universe, while in [11] the use of only Friedmann spacetimes 
leads to a universe with q — ► + at late times. 
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